When Do Like Charges Attract? 
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We study the interaction potential between two polyions inside a colloidal suspension. It is 
shown that at large separation the interaction potential is purely repulsive, with the induced attrac- 
tive interactions being doubly screened. For short separations the condensed counterions become 
correlated, what leads to an effective attraction between the two macromolecules. 
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The systems in which the interactions between the par- 
ticles are predominantly due to the long-ranged Coulomb 
force remain an outstanding challenge to classical statisti- 
cal mechanics. Even such basic question as the existence 
of a phase transition in a symmetric electrolyte has re- 
mained uncertain until quite recently jjj. Although the 
answer to this question has proven to be affirmative Q], 
the true nature of the transition as well as its universal- 
ity class still remain unclear |3|. Surprisingly the theory 
that is in closest agreement with the Monte Carlo simula- 
tions is based on the old ideas first introduced by Debye, 
Hiickel, and Bjerrum Q more than seventy years ago. 
The fundamental insight of Debye was to realize that 
since the mean force inside the electrolyte is zero, it is 
the correlations in positions of oppositely charged ions 
that produce the main contribution to the free energy. A 
great advantage of the Debye-Huckel (DH) theory, be- 
sides its simplicity, is that its linear structure allows it to 
avoid the internal inconsistencies that are often present in 
the more complicated non-linear theories of electrolytes 
||. The price for linearity, however, is that the DH 
theory cannot account for the non- linear configurations, 
such as formation of dipoles, which become important 
at low temperatures. It was an idea of Bjerrum, pro- 
posed only three years after the publication of original 
DH paper, that the missing non-linearities can be rein- 
troduced into the DH theory through the assumption of 
chemical equilibrium between the monopoles and mul- 
tipoles. The extended Debye-Huckel-Bjerrum (DHBj) 
theory has proven extremely successful in elucidating the 
underlying physics of symmetric electrolytes ||, polyam- 
pholytes || , rod- like polyelectrolytes , and charged col- 
loidal suspensions Q , with its validity extending far into 
the regime where the pure linearized DH theory fails. In 
this letter we shall use the DHBj theory to explore one 
of the most fascinating phenomena in condensed matter 
physics, the appearance of attraction between two like- 
charged polyions inside a colloidal suspension ||. 

To model a charged colloidal suspension we use a 
restricted primitive model, in which the polyions are 
treated as hard spheres of radius a and uniform sur- 
face charge, tr_ = — Zq/Aita 2 , and the counterions are 
point particles of charge zq. The solvent is modeled as a 
uniform medium of dielectric constant D. In Ref. B it 



was demonstrated that the equilibrium state of the col- 
loidal suspension consists of some free polyions of den- 
sity poj free counterions of density pf, and of clusters 
composed of one polyion and of 1 < n < Z/z associ- 
ated counterions. The effective charge of a n-cluster is 
Z eS = -(Z - nz)q = 4na 2 Ao-. The DHBj theory al- 
lows us to explicitly calculate the distribution of cluster 
densities {p n } M. 
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FIG. 1. The two clusters separated by a distance 
R — 2a + L. The polyions have a fixed surface charge density 
a- and the mobile surface charge a+ due to the condensed 
counterions, +. The overall imbalance of charge, |<r_| > |cr + 
is responsible for the polarization of clusters and for decrease 
in the electrostatic energy relative to the uniform distribution 
of charge, Act — <j+ + ct_ . 

Consider two clusters separated by a distance R = 
2a + L, each consisting of a polyion and n condensed 
counterions. It is convenient to set up a bipolar coor- 
dinate system with the two polyions located at r% = Q 
and r 2 — 0, respectively. To calculate the effective po- 
tential between the two clusters we shall appeal to the 
DHBj theory The electrostatic potential inside the 
suspension satisfies the Poisson equation, 



V 2 c6 = 



0, 

-47r / 9 g (r 1 ,r 2 )/£), 



for r± < a or r 2 < a, 
for n > a and r 2 > a. 



(1) 



Within the DHBj theory the charge density is approxi- 
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mated by 



Pg(ri,r 2 ) 



Z/z 
n=0 



zn)qp n + zqp f e 



-l3zqtj>(r 1 ,T 2 ) 



+ a-8(n - a) + a+e-P^^^Sin - a) (2) 
+ cr_5(r 2 - a) + a+e-P zq ^ ri ^5{r 2 - a), 

where, (3 = 1/fcsT, and the average surface charge den- 
sity of condensed counterions is c + = nzq/Aira 2 . Note 
that only the counterions are assumed to get polarized, 
while the polyions and clusters contribute only to the 
neutralizing background. As a second approximation, 
linearization leads to p q {rx,r-i) = — (Z)/47r)K 2 0(ri, r 2 ) + 
[Act — Da0(ri, r 2 )/27r][(5(ri — a)+(5(r 2 — a)], where the in- 
verse Debye and the Gouy-Chapman lengths are respec- 
tively k = y/ '4-7T (3 'z 2 [ q 2 p j , and a = 2ito + f3zq/ D. Naively 
one might think that linearization is valid only at high 
temperatures. This, however, is not the case, since the 
non-linearities are effectively included in the renormal- 
ization of the polyion charge by the formation of clusters 
|8|,|l0|| . The set of equations (|l|) is extremely difficult to 
study due to non-trivial boundary conditions. Neverthe- 
less some progress can be made in the two limiting cases, 
L ^> a and L -c a. All the detail of calculations will 
be presented elsewhere, here we shall content ourselves 
with giving a simple and an intuitive explanation of the 
results. 

In the limit L»a and kL > 1 one finds that the cru- 
cial small parameter is e = exp(—nR)/R |11|. It is then 
possible to show that the interaction potential between 
the two n-clusters is 



w(R) = z 2 s e 2 (Kay 



DR 



(na)h(aa)- 



-2kR 



DR 2 



(3) 



where the two scaling functions are 9{x) — e x /(l + x) 
and h(x) = | — j- ln(l + The first term is just 

the usual DLVO potential between the two colloidal par- 
ticles with the effective charge Z e g The second 
term is the result of two types of interactions. First, 
presence of polyions produces holes in the ionic atmo- 
sphere. Since the charges of the holes have the same sign 
as the charges of the polyions, the holes interact repul- 
sively with the polyions JTlJ . The second effect is that 



the mobile counterions confined to the surface are easily 
polarized. Thus, the electric field produced by the cluster 
C\ induces a dipole moment in the cluster C 2 , and vice- 
versa. The polarizability of a n-cluster is found to be 
7 = 2aa 3 /(3 + lad). The effective dipole moment inter- 
acts attractively with the charge that induced it. Since 
the scaling function h{x) changes sign, we see that the 
correction to the DLVO potential is repulsive (hole dom- 
inated) for aa < 1.716... and is attractive (dipole domi- 
nated) for aa > 1.716.... However, we do note that this 
induced attraction is "doubly" screened and is dominated 
by the leading order DLVO term. 

We shall now turn our attention to the opposite limit 
L«n and kL < 1. Under these conditions the relevant 
small parameter is e = L/a, and the Debye screening due 
to unassociated counterions can be neglected. It is con- 
venient to define a cylindrical coordinate system, (p,y), 
with the y axis passing through the centers of the two 
clusters, Fig 1. The two polyions are located at y = —a 
and y = a + L, respectively. The interaction potential 
can be subdivided into two parts. The mean-field repul- 
sion, VmFi arising from the overall imbalance of charge 
on the two clusters, and the attraction coming from the 
correlations among the condensed counterions. It is evi- 
dent that the repulsive energy satisfies, Vmf < Z 2 S /DR. 
In what follows we shall approximate Vmf by the upper 
bound, Vmf = Z 2 ff/DR, that is, we shall treat the mo- 
bile charge, cr+, as, on average, uniformly distributed on 
the surface of the two polyions. Clearly this assumption 
will overestimate the importance of repulsion, since it 
neglects the ability of the bound counterions to arrange 
themselves in the most efficient way to minimize the en- 
ergy, see Fig. 1. With this in mind, to leading order in 
e, the geometry can be replaced by that of two plates 
located at y = and y = L, each with a fixed surface 
charge, and the mobile charge, The free energy 
can be calculated using the DH theory. Let us fix one 
counterion at (0,0). The equations (|l|) may now be in- 
tegrated to give 



,. . 7rAcr , 2irAa ,. . . , . . 

4>{P, V) = —f^L —{\y\ + \ y -L\) + A</> c (p, y), 



D 



D 



(4) 



where the correlation potential can be written in terms 
of the zero-order Bessel function, 



J °° A(k)e ky J (kp)dk, fory<0, 
Ac/> c (p,y) = \ Jv°[B(k)e k y +C(k)e- k y]J Q (kp)dk, for < y < L, 

J °° E(k)e- k yj (kp)dk, for y > L. 



(5) 
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The coefficients, A, B, C, E, are determined through the 
conditions of continuity of the potential and disconti- 
nuity of the normal component of the electric field re- 
lated to the presence of surface charge at y = and 
y = L, u(p,0) = Act — a + zq/3(p(p 1 0) + zq5(p)/2irp 
and cr(p,L) = Act — a + zq/3(f>(p, L), respectively. For 
our purpose it is sufficient to determine only the co- 
efficient A(k), since this is the only one that will en- 



ter into the calculation of free energy. We find A(k) = 
[k(k + a) - kae- 2kL ]/[(k + a) 2 - a 2 e - 2kL }. The electro- 
static potential of the counterion due to the presence of 
other charges is lim p ^ ((j)(p, 0) — zq/Dp) fl3]j . The regu- 
larized free energy per unit area is then obtained through 
the Debye charging process, in which all the particles in 
the system are charged from to their final charge, 



2n(Aa) 2 
fn = £) 



I MX+ D + I MX l _{k + X 2 a) 2 -X i a 2 e~ 2kL k + X 2 , 



o 



kdk 



(6) 



The force per unit area that each polyion exerts on the 
other is F = —dfn/dL. Performing the differentiation 
wc find F = 2l^L> — 1_ ^fe§^, with the scaling function 
g{x) plotted in Fig 2. The first term is the repulsion 
due to excess charge on the two clusters, while the sec- 
ond term is the correlation-induced attraction |fL4| . For 
x ^> 1, 3(00) = — C(3)/87T, where £ is the Riemann zeta 
function. For small x, g(x) rj — x 2 /4ir, and we note that 
for short enough separation the correlation-induced at- 
traction will always dominate the excess charge repul- 
sion! Furthermore, as was mentioned earlier, taking a 
better account of the spherical geometry will only favor 
the attraction by decreasing the repulsive contribution, 
Vmf-i an d increasing the correlations (a) between the in- 
duced counterions on the two hemispheres facing each 
other, Fig 1. Let us now make some estimates. Suppose 
that the condensed counterions neutralize a fraction of 
fixed charge, |ct+| = We shall see attraction when 

/i = z 2 Ab/ 2 /(1 — f) 2 > L, where the Bjerrum length, 
Xb = fiq 2 /D, is 7.2 A in water. Consider a colloidal sus- 
pension of polystyrene particles with characteristic size 
700 A and charge Z = 1000. In Ref. @ it was found 
that the average cluster size is < n >= 400 or / = 0.4, 
what leads to p = 3.2 A. This is smaller than the size 
of a hydrated counterion. We, therefore, do not expect 
that for colloids with Z < 1000 one will observe any at- 
traction. However, for polyions with Z = 3000, / = 0.73, 
the attractive interaction will dominate for L < 53 A and 
may be observable experimentally. If a multivalent salt 
is added to the colloidal suspension, due to a strong elec- 
trostatic attraction the multivalent counterions will be 
preferentially adsorbed to the polyion surface, leading to 
even stronger fluctuation-induced attraction. We like to 
emphasize that a presence of attraction does not imply 
the existance of a phase transition. The thermodynamic 
properties of a colloidal suspension are mostly deter- 
mined by the counterions and their interaction with the 
polyions and clusters. The contribution of the cluster- 
cluster interaction to the osmotic pressure is minimal 
M. Nevertheless the metastable effects connected with 
the presence of attractive forces might explain the un- 



usual observations connected with the charged colloidal 
suspensions Finally, we should stress that the at- 
traction predicted by the DHBj theory is intrinsically a 
finite-concentration result and will disappear at infinite 
dilution. The reason for this is that, in the case of spher- 
ical colloids, the formation of clusters is a purely finite- 
concentration phenomena. This should be contrasted 
with the rod-like polyelectrolytes, for which the strong 
logarithmic potential existing between the polyions and 
the counterions allows for the condensed layer to persist 
all the way to zero density [Q, making it viable to study 
just two polyions with their counterions Jig] . 
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FIG. 2. The scaling function g(x). 
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